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Abstract 

We investigate the possibility that the broken spatial inversion symmetry by a trap potential 
induces a spin-triplet Cooper-pair amplitude in an s-wave superfluid Fermi gas. Being based 
on symmetry considerations, we clarify that this phenomenon may occur, when a spin rotation 
symmetry of the system is also broken. We also numerically confirm that a triplet pair amplitude is 
really induced under this condition, using a simple model. Our results imply that this phenomenon 
is already present in a trapped s-wave superfluid Fermi gas with spin imbalance. As an interesting 
application of this phenomenon, we point out that one may produce a p-wave superfluid Fermi 
gas, by suddenly changing the s-wave pairing interaction to a p-wave one by using the Feshbach 
resonance technique. Since a Cooper pair is usually classified into the spin-singlet (and even-parity) 
state and the spin-triplet (and odd-parity) state, our results would be useful in considering how 
to mix them with each other in a superfluid Fermi gas. Such admixture has recently attracted 
much attention in the field of non-centrosymmetric superconductivity, so that our results would 
also contribute to the further development of this research field, on the viewpoint of cold Fermi 
gas physics. 

PACS numbers: 03.75.Ss, 03.75.-b, 67.85.Lm 


1 


I. INTRODUCTION 


Since any experiment on a superfluid Fermi gas is done in trap potential[l-7|, it is in¬ 
inhomogeneity. An 


Another example 


teresting to explore physical phenomena originating from this spatia’ 
example is surface oscillations observed in a ®Li superfluid Fermi gasjs, 
is the phase separation observed in a ®Li Fermi gas with spin imbalance [5l-l7|, where the 
spin-balanced superfluid region in the trap center is spatially surrounded by excess atoms. 

In addition to these macroscopic phenomena, the spatial inhomogeneity can also affect 
microscopic superfluid properties. Noting that a trap potential breaks the spatial inversion 
symmetry when the inversion center is taken to be away from the trap center, we expect that 
the parity becomes no longer a good quantity to classify the spatial structure of a Cooper 
pair, leading to the admixture of even and odd parity symmetry. Since a pair wavefunction 
is always antisymmetric with respect to the exchange of two fermions, this naturally leads 
to the mixing of spin-singlet and spin-triplet state. When this phenomenon occurs, the 
s-wave superfluid state is accompanied by a triplet Cooper pair amplitude, in addition 
to the ordinary singlet component. (The Cooper pair amplitude is symbolically written as 
{cp^aC-p,a'), where Cp^a is an annihilation operator of a Fermi atom with pseudospin a =t, i.) 

The purpose of this paper is to theoretically explore this possibility in a trapped s-wave 
superfluid Fermi gas. Using symmetry considerations, we prove that this phenomenon may 
occur, when a spin rotation symmetry of this system is also broken, in addition to the bro¬ 
ken inversion symmetry by a trap potential. In a two-component Fermi gas, this additional 
condition is realized, when two species feel different trap potentials or chemical potentials, 
or when they have different atomic masses. Although this is a necessary condition, we nu¬ 
merically conhrm that a triplet pair amplitude is really induced under this condition, within 
the mean-held theory for a model two-dimensional lattice Fermi superhuid in a harmonic 
trap. 

In considering a triplet pair amplitude, one should note that the appearance of this 
quantity does not immediately mean the realization of a triplet superhuid state. Actually, 
the system is still in the s-wave superhuid state, as far as the system only has an s-wave 
interaction. This is simply because the symmetry of a Fermi superhuid is fully determined 
by the symmetry of the superhuid order parameter, which is essentially given by the product 
of a pairing interaction and a pair amplitude. For example, an s-wave superhuid Fermi gas 
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with a contact type s-wave pairing interaction —Ug (< 0) is characterized by the ordinary 
s-wave superfluid order parameter, 


Us ^ y ( 1 ) 

P 

which is flnite when the pair amplitude (cp^^c-p^i) has the s-wave component. The odd-parity 
component does not contribute to in Eq. ([T]). 

However, for an s-wave superfluid Fermi gas with both the singlet and triplet pair ampli¬ 
tude, when one suddenly changes the s-wave pairing interaction to a triplet (and odd parity) 
one U{p,p'), while the s-wave superfluid order parameter in Eq. ([1]) immediately vanishes 
due to the vanishing s-wave interaction {Ug = 0), the product of the triplet interaction and 
the triplet component in the pair amplitude {cp^^c-p^i) (which is assumed to have already 
existed in the s-wave state) immediately gives a flnite triplet superfluid order parameter. 


A(p) = ^f/(p,p')(cp',tc-p',t)> (2) 

p' 

when the triplet interaction U{p,p') is chosen so that the momentum summation in Eq. ([2]) 
can be flnite. In an ultracold Fermi gas, the change of the interaction is possible by using 
a tunable interaction associated with a Feshbach resonance js 


15| . Then, by definition, one 


obtains a triplet superfluid Fermi gas characterized by the superfluid order parameter A{p) 
in Eq. ([2]), at least just after this manipulation. This makes us expect that, when one can 
induce a p-wave pair amplitude in an s-wave superfluid Fermi gas, a p-wave superfluid Fermi 
gas may be realized. This possibility has recently been discussed by one of the ant hors |l6|. 


where a p-wave pair amplitude is induced by a synthetic spin-orbit interaction 
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25j. The 


present paper provides another source of p-wave pair amplitude, without using an artificial 
gauge held. 

The admixture of singlet and triplet Cooper pairs has recently attracted much attention 
in the held of non-centrosymmetric superconductivity 26l. l27l |. where a crystal lattice with 
no inversion center causes this phenomenon. In this held, it has been pointed out that this 
admixture may be the origin of the anomalous temperature dependence of the penetration 


depth observed in Li 2 Pt 3 B 


Thus, an s-wave superfluid Fermi gas with a triplet pair 


amplitude would be also helpful to the study of this electron system. 

This paper is organized as follows. In Sec.II, we clarify the necessary condition for a 
triplet Cooper pair amplitude to appear in a trapped s-wave superfluid Fermi gas. In 
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Sec.Ill, we numerically evaluate how large a triplet pair amplitude is induced under the 
condition obtained in Sec.11. In this section, we treat a superfluid Fermi gas loaded on a 
two-dimensional square lattice, within the mean-held theory. Throughout this paper, we 
take h = /cb = 1, for simplicity. 


II. CONDITION FOR TRIPLET PAIR AMPLITUDE TO APPEAR IN A 
TRAPPED s-WAVE SUPERFLUID FERMI GAS 


We consider a three-dimensional s-wave superhuid Fermi gas, described by the Hamilto¬ 
nian, 


H 





f/,T|(r)T|(r)T^(r)Tt(r) 


( 3 ) 


Here, TQ,(r) is a fermion held operator with pseudospin a =t, describing two atomic 
hyperhne states. —Ug (< 0) is a contact-type s-wave pairing interaction. ha^a'{r) is a one- 
particle Hamiltonian density, consisting of a kinetic term and a potential term, detailed 
expression of which will be given later. 

We assume that the system is in the ordinary s-wave superhuid state with the s-wave 
superhuid order parameter. 


A,(r) = f/,(T^(r)T|(r)). 


( 4 ) 


We also assume that any other spontaneous symmetry breaking is absent (such as the triplet 
superhuid state). 

In this model superhuid, we consider the spin-triplet Cooper-pair amplitude, given by 




(vl/^(r)vl/^(r')) = 1), 

^ [(Tt(r)vl/^(r')) + (vl/^(r)vl/^(r'))] (F. = 0), 

(T^(r)vl/^(r')) = -1), 


( 5 ) 


\ 


where denotes the ^-component of the total spin of each pair amplitude. The triplet pair 
amplitude in Eq. ([5]) does not contribute the s-wave superhuid order parameter A 5 (r) in 
Eq. dH), because <I)f^(r,r) = 0. The spin-singlet pair amplitude. 


$s(r,r') 

only contributes to Eq. fjTl) . 
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[(Tt(r)T^(r')) - (T^(r)Tt(r'))], 


( 6 ) 
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We first prove that the broken spatial inversion symmetry is necessary for a triplet pair 
amplitude to appear in an s-wave superfluid Fermi gas. For this purpose, we conveniently 
introduce the inversion operator P{R) with respect to the inversion center R. The held 
operator is transformed under this operation as 


T«(r) = P(i^)T„(r)p-l(i^) = T«(i^ - /), 
where r = R + 1. The inverted Hamiltonian H = PHP~^ is then given by 


( 7 ) 


( 8 ) 


- U,^{R-l)^l{R-l)^^{R-l)^^{R-l) 

When the one-particle Hamiltonian density ha^a'i^) has the symmetry ha^a'{R + 1) — 
ha,a>{R — 1), this system is invariant {H = H) under this symmetry operation. On the 
other hand, the triplet pair amplitude in Eq. ([S]) with the center of mass posi¬ 

tion R= [r + r']/2 is transformed as, 

l>f=i(r,r') = (Tt(r)Tt(r')) = (^t(-R-W2)^t(^ +W2)) 

= + r,el/2)Tt(i^ - rrei/2)) = -<Ff='(r, r'). 


(9) 


where Vj-^i = r — r' is the relative coordinate. We also hnd <F 


5^=0,-!, 


r,r') = -<F 


S^=0,-1/ 


r, r' 


That is, the triplet pair amplitude <Ff^(r,r') vanishes, when the system has the inversion 
symmetry {H = H) with respect to the inversion center R = [r + P\l2.. Thus, the broken 
inversion symmetry is necessary for a triplet pair amplitude to appear. 

For the singlet pair amplitude in Eq. ([6]), this symmetry operation simply gives 
<hs(r,r') = $s(r,r'). As expected, this quantity may be hnite. 

The one-particle Hamiltonian density hQ, Q,/(r) in the ordinary uniform Fermi gas has the 
form 

r 1 

ha,a' F ^a,a'1 (^0) 

where p = —iV, m is an atomic mass, and /i is the Fermi chemical potential. Equation fITU]) 
has the symmetry property, ha,a'{R + 1) = ha,a'{R — 1), with respect to I for an arbitrary 
R. To conclude, any triplet pair amplitude is not induced. 

In the presence of a harmonic trap, the one-particle Hamiltonian density becomes inho¬ 
mogeneous as 


hn 


y'[r) = 


PL 

.2m 


h + 


( 11 ) 
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so that it does not have the inversion symmetry except at -R = 0. However, when we 
consider the s-wave superfluid state in this trapped case, any triplet pair amplitude is not 
actually induced (although we do not explicitly show the result here). Of course, since 
the condition obtained from the inversion symmetry is a necessary condition, the broken 
inversion symmetry does not guarantee the appearance of a triplet pair amplitude. 

In this regard, we point out that the vanishing triplet pair amplitude in the trapped case 
is due to the fact that this system still has a rotation symmetry in spin space. To see this, 
we next consider the spin rotation of the held operator, given by 


T„(r) = R{e)^^{r)R-\e) = T„,(r). (12) 

V / a,a' 
a' 

Here, 6 = 6eo describes the spin rotation around the unit vector eg with the angle 6, 
and (T = {(Tx,o'y,az), where aj (j = x,y,z) are Pauli matrices. (As usual, we take the 
spin quantization axis in the a^-direction.) For the three “tt rotations” specified hy 0 = 
(7r,0,0) (= 9 ^), (0,7r, 0) (= 0^), (0,0, tt) (= 0^), Eq. (IT^ can be written as. 


I ^t(^) 



= iaj 



(13) 


where we have used the formula = cos(0/2) + iajSm{6/2). Under the vr-rotation, the 
Hamiltonian in Eq. ([3]) is transformed as 


H = R{ei)HR-\ei) 

= f dr 

a 

a ’ 

Here, h{r) = Thus, one has H = H, when 

[h{r),aj] = 0 


(14) 


(15) 


is satisfied. 

While the singlet pair amplitude in Eq. ([6]) remains unchanged under these vr-rotations. 
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the triplet component is transformed as 















if-”(r,r') 

= < 


1 

'll 

0 

(if-—(r.r')) 











< 



(9 = »%), 

(9 = 9;), 

(9 = 9;). 


( 16 ) 


For example, when we set 6 = 6^, Eq. flTHl) means that r') = 0, when [h{r), ax] = 0. 

(The other two components with Sz = ±1 are not excluded in this case.) When the system 
invariant under all the vr-rotations ( 0 ^’^’^), any triplet pair amplitude is not induced. 

To conclude, the broken spin rotation symmetry characterized by is necessary for a 
triplet pair amplitude to be induced in a trapped s-wave superfluid Fermi gas. This is the 
reason why the model case described by Eqs. (El) and flTTll is not accompanied by any triplet 
pair amplitude. 

The broken spin rotation symmetry is realized, when the strength of trap potential K 
in Eq. (ITTD depends on spin (= iFa). In this case, the one-particle Hamiltonian h{r) = 
{ha,a'{r)} can be written as 


h{r) = 


, K^ + K^ 2 


At-A'l 2 

+ - -r^cr,. 


( 17 ) 


Equation (lT5l) is then satished only when j = z. Thus, we hnd from the last line in Eq. (1T6|) 
that only the triplet pair amplitude with = 0 may be induced. Since the last term in Eq. 
(IT7|) works as an external i nag netic field, this phenomenon is also expected in the presence 


of spin imbalance js-I t], [2^, [3^, where two species feel different Fermi chemical potentials 


p-|- 7 ^ fii- Another possibility is a trapped hetero superfluid Fermi gasj3ll-l39|. where two 
species have different atomic masses In Sec. Ill, we will numerically examine 

these cases. 

Before ending this section, we briefly note that the broken inversion symmetry, as well as 
the broken spin rotation symmetry, are also realized in a spin-orbit coupled s-wave superfluid 
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Fermi gasjl7H25j. Indeed, Refs. 
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41| predict that a p-wave pair amplitude is induced 


in this case. Although we do not deal with this case in Sec.Ill, we explain in Appendix A 
how to apply the present symmetry consideration to this case. 


III. NUMERICAL CONFIRMATION FOR THE INDUCTION OF TRIPLET 
PAIR AMPLITUDE IN A TRAPPED s-WAVE SUPERFLUID FERMI GAS 


To examine whether or not a triplet pair amplitude is induced under the condition ob¬ 
tained in Sec. II, we consider a model s-wave Fermi superfluid loaded on a L x L two- 
dimensional square lattice, within the mean-held approximation. Although this simple model 
cannot be directly applied to a real continuum Fermi gas, it is still helpful to grasp basic 
characters of this phenomenon. 

The Hamiltonian is given by 


HmF — — ta 


+ h.C. 




+ h-c- 


E - Usn_a{ri) 




(18) 


2,0 


Here, cj.. „ is a creation operator of a Fermi atom at the lattice site r, = with 

pseudospin a =t, and the Fermi chemical potential /Iq. —ta describes a particle hopping 
between the nearest neighbor sites, and (q j) means the summation over the nearest neighbor 
pairs. In Eq. (ITSH . the s-wave superhuid order parameter As(rj) = f^s(cri,tCr.._ 4 ,), as well 
as the Hartree potential —UsU-airi) = — -a), are obtained from the mean-held 
approximation for the on-site Hubbard interaction —where —Ug (< 0) 
is the interaction strength. Ua(rj) = is a harmonic trap potential, where is 

the strength of a trap potential which a-spin atoms feel. Here, the spatial position r* is 
measured from the center of the L x L square lattice, and is the distance between the 
trap center and the edge of the system. For simplicity, we take the lattice constant a to be 
unity. 

In the present model, the spatial inversion symmetry is broken by the trap potential 
except at the trap center. For the spin rotation symmetry, it is broken when one of ta, 
/Tq, and Vf^, depends on pseudospin a =t, Since all these cases satisfy Eq. (lT5j) only 
when j = z, Eq. ffT6l) indicates that one may only consider the possibility of the triplet pair 















amplitude with Sz = 0. In the present model, this component is given by 


$ 


Sz=0, 


r^rj) = 


^/2 




(19) 


For comparison, we also consider the ordinary singlet component, given by 




V 2 I 


( 20 ) 


Besides the superfluid order parameter, the condensate fraction is also a fundamental 


quantity in the superfluid phase 
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45| . which physically describes the number of Bose- 


condensed Cooper pairs. In an ordinary s-wave superfluid state, it has the form. 


m = 




( 21 ) 


R={ri+rj)/2 

where the local condensate fraction n®(i?) is directly related to the singlet pair amplitude 
in Eq. (|20l) as 

<(^) = ^ E l$s(i^ + W2,i^-rrel/2)|2. (22) 

rTel=ri-rj 

In addition to the singlet component of the condensate fraction iV® in Eq. fl?I|) . the present 
system may also have the spin-triplet component [la] , iV* = X)K={ri+rj)/ 2 where 


1 


y(-R) = ^ E W2)i' 


(23) 


"^rel —j 


The total condensate fraction is given by + N^. In what follows, we simply call iV® and 
Ad the singlet and triplet condensate fraction, respectively. 

We note that the square lattice in our model does not affect the symmetry consideration 
in Sec. II. This is because the square lattice is invariant under the spatial inversion with 
respect to the center of mass position R = [r, -|- rj\/2 of a triplet pair amplitude <Ff"'(rj, Vj). 
In addition, since the spin rotation symmetry is also unaffected by the crystal lattice, the 
necessary condition obtained in Sec. II is still valid for the present case. 

As usual, we diagonalize the mean-field Hamiltonian Hmv in Eq. ffTSj) by the Bogoliubov 


transformation 


46l |. Since this is a standard procedure [47|, we do not explain the detail 


here, but summarize the outline in Appendix B. We numerically carry out the Bogoliubov 
transformation, to self-consistently determine As{ri), Uaivi), and /i„. We then evaluate the 
triplet pair amplitude Tj) in Eq. ([19]). 
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FIG. 1: (Color online) (a) Calculated triplet pair amplitude 0), (rj, + 1,0)) along the x- 

axis {rjy = 0) in an s-wave superfluid Fermi gas with trap-potential difference, (b) s-wave superfluid 
order parameter As(r^,0). (c) Density profile nQ,(r*,0). We take Vq/Vq = 0.5, t^/ti = 1, and 
Us/t = 2.5. This parameter set is also used in Figs. [2]l4l 


In numerical calculations, we take the lattice size L = 41, and Vq = 2t, where t = + 

ti\/2. To avoid lattice effects, we consider the low density region, by setting = 59 

in the absence of spin imbalance (where is the number of Fermi atoms in the a-spin 

component). The total number N = of Fermi atoms then equals N = 118. In this 

case, the particle density is at most Uaivi) < 0.3 1 even in the trap center. We take a 

small but hnite temperature T/t = 0.01, in order to suppress effects of discrete energy levels 
associated with the hnite system size. 
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FIG. 2: (Color online) Calculated triplet pair amplitude r^) in an singlet superfluid Fermi 

gas with trap-potential difference, (a) (r* -|- l,rp). (b) , r^, (r^,r* -|- 1)). 

(c) {rl + l,ry + 1)). (d) Singlet pair amplitude <l>s(rj,ri). 


Figure [11(a) shows the evidence that the triplet pair amplitude with 5*2 = 0 is induced in 
the s-wave superfluid state, when both the spatial inversion symmetry and the spin rotation 
symmetry are broken by the trap potential Va(rj). From the comparison of this hgure with 
FigslUb) and (c), one hnds that $f"’^°((r^ + l, 0 ),(r^, 0 )) appears everywhere in the gas 
cloud where the s-wave superfluid order parameter As(r^, 0 ), as well as the atom density 
nQ,(rJ,, 0), are hnite, except at the trap center. Since the system still has the spatial inversion 
symmetry at the trap center, the node structure seen in Fig. [11(a) agrees with the symmetry 
consideration in Sec. II. We emphasize that the triplet pair amplitude is not induced when 
Vq = I/q*', although we do not explicitly show the result here. 

Figure [21(a) shows that the point node seen in Fig[Tl(a) is actually a line node along the 
2 /-axis. This node structure comes from the symmetry property that, while the present 
square-lattice model has the reflection symmetry with respect to the y axis, the triplet pair 
amplitude <Ff"^(rj,rj) behaves as, 

+ r,ei/2, R - r,ei/2) = -$^(i^ - Cei/2, R + (24) 

when R= [ri + rj\/2 = (0, Ry) and r^ei = — Vj = (r^gj, 0). Since the present lattice model 

is also invariant under the reflections with respect to the x axis, as well as the lines along 
y = ±x, the triplet pair amplitude <I)f"'(rj, r^), with the relative vector rrei = ri — Vj being 
perpendicular to one of them, has the line node along the reflection line. (See Figs. [2](b) and 
(c).) In a continuum system with no lattice, the triplet pair amplitude is expected to always 
have the line node, which is perpendicular to the relative vector of the pair amplitude. We 
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FIG. 3: (Color online) (a) Triplet pair amplitude as a function of the relative coor¬ 

dinate Trei = Vi — Tj. We take /Vq = 0.5. For each center of mass position R = [ri + rj]/2 
(open square), the pair amplitude is plotted inside the region between two dashed line, by taking 
R as the origin, (b) Spatial variation of synthesized p-wave superfluid order parameter Ap{p, R). 
In this panel, the p-dependence of Ap(p, R) is schematically shown, being centered at R (open 
square). 


briefly note that such a node is not obtained in the singlet component, as shown in Fig. 

Eld). 

Figure [3la) shows the spatial structure of the triplet pair amplitude r^) with 

respect to the relative coordinate Urei = r, — Vj. Noting that the pairing symmetry is 
determined by the angular dependence in relative-momentum space, we hnd that the induced 
pair amplitude has the p-wave symmetry. That is, the pair amplitude has the Pa;-wave {py- 
wave) symmetry, when the center of mass position is on the x-axis (p-axis). 

An advantage of the cold Fermi gas system is that one can tune the pairing interaction by 
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adjusting the threshold energy of a Feshbach resonance. Although this technique is usually 
used to adjust the interaction strength for a fixed interaction channel, one may also use 
this technique to change the interaction channel from the s-wave one to a p-wave one. In 


an s-wave superfluid Fermi gas with triplet pair amplitude s’ 
suddenly change the s-wave interaction to the p-wave one 48 


lown in Fig. Ela), when one 


49l, 


hfp-wave — Up p ■ p 

P,P',Q 


(25) 


the p-wave superfluid order parameter, 


Ap{p,R) = f/p5^p■p'<Ff=0(p',i^), (26) 

P' 

would immediately become hnite. Here, R) is the Fourier-transformed triplet pair 

amplitude with respect to the relative coordinate rrei. We emphasize that this triplet pair 
amplitude has already existed before the change of the interaction. Thus, just after this ma¬ 
nipulation, we expect the spatial structure of the induced p-wave superfluid order parameter 
schematically shown in Fig. |3](b). The s-wave superfluid order parameter immediately dis¬ 
appears because of the vanishing s-wave interaction {Us = 0), and the s-wave pair amplitude 
‘hs(v’i, Tj) only remains. Thus, at least just after this manipulation, by definition, the system 
is in the p-wave superfluid state with the synthesized p-wave superfluid order parameter in 
Eq. (I26il . This unconventional superfluid phase would be in the non-equilibrium state, so 
that we need further analyses on the time evolution of this state. However, the combined 
Feshbach technique with the induced triplet pair amplitude is an interesting idea to realize 
a p-wave superfluid Fermi gas. 

Figure m shows the local condensate fraction rf^^{R) in an s-wave superfluid Fermi gas 
with trap-potential imbalance. In panel (a), the triplet component nl{R) is enhanced around 
\R\ = 6, as well as the region near the trap center (except at i? = 0, where the triplet con¬ 
densate fraction vanishes). On the other hand, panel (b) shows that the singlet component 
nl{R) monotonically decreases, as one goes away from the trap center. The latter behav¬ 
ior is consistent with the spatial variation of the s-wave superfluid order parameter As(rj) 
shown in Fig. [11(b). 

The large triplet condensate fraction nj;(i?) near the trap center seen in Fig. |4](a) is due to 
the spin imbalance (n'f(rj,, 0) > ni^rl, 0)) in the trap center. (See Fig. [11(c).) This naturally 
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FIG. 4: (Color online) Local condensate fraction in an s-wave superfluid Fermi gas with trap- 
potential difference, (a) Triplet component The intensity is magnihed to ten times, (b) 

Singlet component 


leads to the broken spin rotation symmetry through the Fermi chemical potential as 
well as the Hartree potential —Usn_a{ri) in Eq. ffTS]) . Thus, although two spin components 
feel almost the same trap potential ( 14 (^«) — around the trap center, the triplet 

condensate fraction is enhanced there (except at R = 0). 

While the difference 14(^*) ~ becomes remarkable as one goes away from the 

trap center, the spin imbalance (n'|-(r*,0) — 'n|(r*,0)) becomes small to eventually vanish 
at \R\ ~ 9. (See Fig. [11(c).) In the outer region, the spin imbalance again occurs as 
n-|-(r^,0) < n|(r40). These enhance ~ 6), as shown in Fig. 111(a). 

Summing up the local condensate fraction nl’^{R) in the gas cloud, one obtains the 
condensate fraction iV®’* in Fig. [5l As expected, panel (a) shows that the triplet component 
Ah is enhanced when Vq /Vq -C 1. We also hnd that Ah becomes large in the intermediate 
coupling regime but becomes small when S> 1. In the strong coupling regime, most 
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FIG. 5: (Color online) Condensate fraction Nc’^ = ™ ™ s-wave superfluid Fermi gas 

with trap-potential difference, (a) Triplet component N^. (b) Singlet component N^. (c) and 
N^, as functions of the interaction strength Ug, when Vq/V^ = 0.5. In panels (a) and (b), the 
region above the solid line is in the superfluid state within our numerical accuracy. (Note that we 
take Tjt = O.OI > 0 in our numerical calculations.) The region between the solid line and the 
dashed line is in the FFLO phase, being characterized by a spatially oscillating superfluid order 
parameter As{ri). These lines are also drawn in FigjGj 

Fermi atoms form singlet molecules, which suppresses effects of broken inversion and spin 
rotation symmetry. Indeed, Fig. 151(b) shows that the singlet component iV® monotonically 
increases with increasing the interaction strength Ug- To clearly see the difference between 
and iV®, Fig. |5](c) shows these quantities as functions of the interaction strength Ug. 
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FIG. 6 : (Color online) Condensate fraction Nl'^ in a trapped s-wave superfluid Fermi gas in the 
presence of mass imbalance t^/t^ / 1. (a) Triplet component N^. (b) Singlet component N^. (c) 
Nl and as functions of the interaction strength Ug/t, when = 0.5. 


In Figs. [5](a) and (b), one sees the FFLO (Fulde-Ferrell-Larkin-Ovchinnikov) phase 50 


53[ |. In this regard, since we are dealing with a two-dimensional lattice model within the 


simple mean-held theory, it is unclear whether or not the FFLO phase still remains in a 


realistic three-dimensional continuum Fermi gas 


54l |. However, since Fig. [5](a) indicates that 


the triplet condensate fraction is also induced in the ordinary BCS region, we hnd that the 
FFLO state is not necessary for the triplet pair amplitude to appear. 

Figure [6] conhrms that the triplet pair amplitude r^) is also induced, when 

the spin rotation symmetry is broken by mass imbalance 7^ !)• In addition. Fig. [7] 
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FIG. 7: (Color online) Condensate fraction A^c’* in a trapped s-wave superfluid Fermi gas with spin 
imbalance / 1. (a) Triplet component N^. (b) Singlet component N^. (c) Ad and N^, as 

functions of the interaction strength Ug/t, when [A^j- — A^j_]/A^ = 0.17 (A^j- — A^^ = 20). In panels (a) 
and (b), the region above the solid line is in the superfluid phase. In the superfluid region shown 
in this figure, the superfluid order parameter in the outer region of the gas cloud always exhibits 
a FFLO-type oscillation in the radial direction. 


shows that this phenomenon also occurs in a trapped s-wave superfluid Fermi gas with spin 
imbalance (A^-|-/A^ 4 , ^ 1) In the latter case, the Fermi chemical potential /i^ depends on 
pseudospin a =t, which breaks the spin rotation svmmetrvj56l|. 

In the presence of spin imbalance, the phase separation is known to occurjs-?!, where 
the superfluid region in the trap center is surrounded by excess atoms. Figure |8](a) shows 
this case. In this panel, since the spin imbalance is almost absent around the trap center. 
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FIG. 8: (Color online) Calculated density profile nQ(rJ,,0) and the triplet 
raj;(i?2:,0) in a trapped s-wave superfluid Fermi gas with spin imbalance. 
Us/t = 3. We set 


condensate fraction 
(a) Us/t = 6. (b) 


the triplet condensate fraction is suppressed there, compared to the case when the phase 
separation does not occur (panel (b)). In addition, the region around the edge of the gas 
cloud is highly spin polarized, so that the triplet pair amplitude is also suppressed there. 
As a result, when the phase separation occurs, the triplet pair amplitude is localized around 
the edge of the gas cloud of the minority component [a =4-), as shown in Fig. [8](b). 


IV. SUMMARY 

To summarize, we have discussed the possibility of inducing a triplet pair amplitude in a 
trapped s-wave superfluid Fermi gas. Using symmetry considerations, we clarified that both 
the broken spatial inversion symmetry and a broken spin rotation symmetry are necessary 
for this phenomenon to occur. We numerically confirmed that a triplet pair amplitude is 
induced when this condition is satished, within the mean-held theory for a two-dimensional 
lattice model. In this conhrmation, we considered the three cases with (1) trap-potential 
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difference, (2) mass imbalance, and (3) spin imbalance. In the first case, we showed that 
the induced triplet pair amplitude is dominated by a p-wave symmetry. Among the above 


three cases, a trapped s-wave superfluid Fermi gas with spin imbalance has been realized 


a 




Thus, our results imply that a triplet pair amplitude is already present in this system, 
although there is no experimental evidence yet. 

Since the symmetry of a Fermi superfluid is fully determined by the symmetry of the 
superfluid order parameter, the induction of a triplet pair amplitude does not immediately 
mean the realization of a triplet Fermi superfluid. In the present case, the system is still in 
the s-wave superfluid state which is characterized by the s-wave superfluid order parameter, 
even in the presence of a triplet pair amplitude. Under this situation, however, when one 
suddenly changes the s-wave pairing interaction to an appropriate p-wave one, the product 
of the p-wave interaction and the triplet pair amplitude that has been induced before this 
manipulation may immediately give a hnite p-wave superfluid order parameter. Since the 
s-wave superfluid order parameter vanishes, by dehnition, we have a p-wave superfluid state, 
characterized by this p-wave superfluid parameter. The change of the interaction would be 
possible by using the Feshbach resonance technique. 

In cold Fermi gas physics, although the realization of a p-wave superfluid state is a 


crucial challenge, current experiments are 
wave interaction, such as three particle loss 


acin^ various difficulties originating from a p- 
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59l | , as well as the dipolar relaxation ^ . In 


this regard, the above idea may avoid these difficulties to some extent, because the triplet 
pair amplitude is prepared in an s-wave superfluid Fermi gas with no p-wave interaction. 
In addition, since we can start from a hnite value of the p-wave superhuid order parameter, 
the system would be in the p-wave superhuid state for a while, until it is strongly damaged 
by the particle loss and dipolar relaxation after the p-wave interaction is introduced. In this 
sense, the induction of a triplet pair amplitude discussed in this paper is important, not only 
as a fundamental physical phenomenon, but also on the viewpoint of the challenge toward 
the realization of a p-wave superhuid Fermi gas. 

In this paper, we have treated a lattice model to simply conhrm the induction of a triplet 
pair amplitude. To quantitatively evaluate this quantity, we need to extend the present 
analyses to a realistic continuum Fermi superhuid. To assess the idea that one produces a 
p-wave superhuid Fermi gas from the induced triplet pair amplitude, it is also important 
to clarify the time evolution of the p-wave superhuid order parameter after the s-wave 
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interaction is replaced by a p-wave one. These problems remain as our future problems. 
Since the pair amplitude always exists in a Fermi superfluid, our results would be useful for 
the study of this fundamental quantity in cold Fermi gas physics. 
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Appendix A: Triplet pair amplitude in a spin-orbit coupled uniform s-wave super¬ 
fluid Fermi gas 


We consider a uniform s-wave superfluid Fermi gas with a spin-orbit interaction. The 
model Hamiltonian is given by 


^ ~ + Hso Us ‘^p+(7/2,t^-p-l-q/2q^-p'+9/2,P^P'+g/2,t- 


(Al) 


p,a 


P,P 


Here, c|, ^ is a creation operator of a Fermi atom with the kinetic energy = j?/(2m) — p, 
measured from the chemical potential p. The antisymmetric spin-orbit interaction H^o has 
the form H. 


H = ^ ^ {t)\c t 

-'-'so / ^ ^P,Q;'^S0 

where /iso(p) = {hff'} is assumed as 


(A2) 


p,a,a' 


hso{p) = ><±[PxO-x +Py(^y] + A, 


a.. 


(A3) 


Here, Aj^ and are spin-orbit couplings. 

We assume that the system is in the s-wave superfluid state with the superfluid order 
parameter = UsJ2p{(^pt^-pi)^ AVe also assume that any other spontaneous symmetry 
breaking is absent. 

In momentum space, the spatial inversion P is described as Cp^a = Pcp^aP~^ = C-p^a- 
Under this operation, each term in Eq. flAlIl is invariant, except for the spin-orbit interaction. 
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which is transformed as 


H„ = PH„p-' = Y. = Y 4 .c>^T'(-pK.<.' = -H„. (A4) 

PyCt,a' p,a^a' 

Thus, the spin-orbit interaction Hgo in Eq. (IA2D breaks the inversion symmetry. 

For the spin rotation R{6), the three vr rotations {6 = 6l, j = x,y,z) corresponding to 
Eq. (fT^ are given by 



= lan 


0=ei 



(j = x,y,z). 


(AS) 


When A_l = 0 and Xz ^ 0 (single component spin orbit interaction), Eq. flAip is not 
invariant under the spin rotations R{0^) and R{0y), because the spin-orbit interaction Hgo 
is transformed as 


Hso = R{O:'nHsoR-\0:’n = A. cla(^T'cp,a' = "A, = -Fso. (A6) 

p,a,a' p,a^a' 

Thus, Eq. (lAlD is invariant only under the vr rotation with 0 = 0^. Noting that the triplet 
pair amplitude 


(Cp'|'C_p'|-) 


{Sz = 1), 


(p) = 1 


^ + (Cp^C-p^)] [Sz 0 ), 


(Cp4,c_pj^) 


is transformed under the three tt rotations as 


{Sz = -1), 




if—(p) t 

4f-»(p) 


$ 


S,=-l 


(P) / 


t (P) 
S^=0 


\ 


= 




-$^=^(p) 
-$f^=l(p) j 
<f)f=-l(p) 

-$f^=0(p) 

<f>f=l(p) j 
-$f=l(p) 

$&= 0 (p) 


(0 = Ql). 


{0 = 01 , 


\ 


[ \ 


\p) j 


{0 = 01 ), 


we find that only ^(p) may be induced. Indeed, Ref. Q shows that it is induced in this 
case. 


(A7) 


(A8) 
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When A_l 7^ 0, the spin-orbit interaction Hgo is not invariant under any tt rotations 
Within this analysis, one concludes that all the trip let p air amplitudes in Eq. flA7p 
may be induced. However, within the mean-held theory. Ref. [1^ shows that the component 
with 5*2 = 0 is not induced when Aj^ 7 ^ 0 and A^ = 0. This is because, in this two-component 
case, the mean-held BCS Hamiltonian, 


^BCS = + ifso + A, ^ 


(A9) 


p,a 


is invariant under the momentum dependent tt spin-rotation which is followed by the U{1) 
gauge transformation, given by 


Cpt 

Cpl 


= e *2 X e* 2 P^''^ 


Cpt 

Cpl 


(AlO) 


where p± = {px,Py)/^Jvi + Vy- this case, the triplet pair amplitude with 5*^ = 0 is 
transformed as so that one hnds = 0 , as obtained in 

Ref. 


Appendix B: Diagonalization of the BCS Hamiltonian in Eq. (|18|) 


The mean-held BCS Hamiltonian in Eq. flTSl) can be diagonalized by the Bogoliubov 
transformation in real space, given by 


Cri,t 


7i ^ 

,t 

^ri.l 

= W 

7l2 


7l2+i 



\ 72L2 / 


(Bl) 


Here, VE is a 2L^ x 2L^ orthogonal matrix, which is chosen so that Huf in Eq. flTSl) can be 
diagonalized as 

2L? 


Hmf - ^ -Ej-ytyj, 
i=i 


(B2) 
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where Ej is a Bogoliubov single-particle excitation energy. After the diagonalization, the 
superfluid order parameter As(rj), as well as the number density riairi) = {c\.. ^Cr^,a), are 
evaluated as, respectively 

2L^ 

A,(r,) = UsJ2 (B3) 

2L2 

n^{r.) = J2wl^f(E,), (B4) 

i=i 


2L2 

= (B5) 

i=i 

where f{E) = l/[e^^^ + 1] is the Fermi distribution function. The number of Fermi 
atoms in the a-spin component is given by 

L2 

A« = ^n«(ri). (B6) 

We numerically calculate Eqs. fIBlD . and flB3P - flB6p . to self-consistently determine Ag^Vi), 
na{ri), and jia- The triplet pair amplitude Tj) in Eq. (fT9D . as well as the singlet 

pair amplitude ^s{ri,rj) in Eq. ([20]), are then calculated as, respectively. 


2L2 


$ 


5^=0/ 


= ^\Wi,kWj+L^kf{-Ek) + Wi+L^^kWj^kf{Ek) 


k=l 


(B7) 


2L2 








(B8) 


k=l 
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